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Abstract 



Using a finite size scaling form for reunion probability, we show numeri- 
cally the existence of a binding-unbinding transition for Directed polymers 
with random interaction. The cases studied are (Al) two chains in 1+1 di- 
mensions, (A2) two chains in 2+1 dimensions and (B) three chains in 1+1 
dimensions. A similar finite size scaling form for fluctuation establishes a 
disorder induced transition with identical exponents for cases A2 and B. The 
length scale exponents in all the three cases are in agreement with previous 
exact renormalization group results. 
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Disorder is often found to give rise to complicated but rich phenomena in nature, an 
excellent example of which is the existence of the spin glass transition and the spin glass 
phase |]J . Needless to say that one has to look beyond the reality, and study simpler models 
that retain only the essential features believed to be responsible for such new events [|TJ. 
Directed polymersQ (DP) in random media or with random interactions |5]|| provide 
a fruitful basis in this regard because of its inherent simplicity. By virtue of mappings to 
nonlinear, noisy surface growth equations [0,0, and other applications, DPs have become 
relevant in a broader context. A crucial result established for directed polymers in random 
media [^] is the existence of a strong disorder "spin glass" type phase at d — 1, and a lower 
critical dimension, d = 2, beyond which a disorder induced transition exists. For the random 
interaction case there is also a weak to strong disorder transition Directed polymers 

are also simpler than the more complex problem of undirected polymers in random medium 
HH. Similarly, many heteropolymer (undirected) problems, especially biopolymers, require 
random interaction [IU|, and DPs with random interaction quite often serve as a simpler 



solvable limit 111 



In this paper, we concentrate on the disorder induced transition in directed polymers 
with random interaction (RANI), to explore the scaling behaviour of fluctuation due to 
quenched disorder. The essential feature here is the mutual short range interaction which 
is random, as might arise when there is a random charge distribution along the length of 
the polymer. A different version of our model appears in the context of wetting phenomena 
in the presence of disorder at d = 1 (T^j. For the RANI model, unlike the random medium 
problem, an exact renormalization group (RG) can be implemented under certain conditions 
IIU . We like to verify the length scale exponents obtained from the perturbative RG. This 
verification is essential for further applicability of RG for the disorder problem to ensure 
that there is no nonperturbative effects. For comparison, one might point out the situation 



for the random medium problem. For the weak to strong disorder transition, RG [|T3| . p~4 
predicts a length scale exponent 2/(d — 2) for d > 2 while numerical calculations [|T5|,|nj for 
d = 3 gave numbers as different as 4.2 and 6.7. We establish here that, unlike the random 
medium case, the RANI model is better controlled. 

For generality let us start with a system of m polymers each of length N with a random 
m body (mth order multicritical) interaction. The model in a continuum formulation is 
given by 




^=oE/ HP +/ dzv m (l + b(z))H5(v u ^), (1) 



i=2 



where rj(z) denotes the d dimensional spatial coordinate of the ith polymer at contour length 
z. The first term denotes the elastic energy part of the Gaussian chains and the second term 
denotes the mutual random contact interaction among the chains. Note that the interaction 
is always at equal length with a random part b(z) in the coupling constant. The randomness 
is dependent only on the length z at which the interaction occurs. This is what one expects 



Arf+l dimensional directed polymer is a random walk stretched in a particular direction with 
fluctuations in the transverse d dimensional space. 



2 



for random monomer distribution along the backbone. More general randomness can arise 



17| but will not be considered here. In the following analysis we shall consider a binary 
distribution for the disorder with equal probability. 

We study this system (actually, a discrete version of this) numerically by generating the 
partition function recursively along the length of the chain. The cases considered are (Al) 
the two chain problem (m = 2) for d — 1, (A2) two chains in d = 2, and (B) three chains 
with three body interaction (m = 3) for d — 1. The reason for choosing the last two cases, 
as explained below, is that d — 2 (d — 1) is the marginal case for the two body (three body) 
interaction problem. Instead of analyzing large systems as has been done in Ref. |7J for 
m = 2 and d — 1, we use finite size scaling to calculate the exponents for various cases. 

The pure counterpart of the disordered model, Eq. [I] is at present well understood 
through the renormalization group analysis [ |T8| , |T9f or by solving the Schrodinger type equa- 
tions for the two body problem p(| . An exact RG shows that there is a bound state for any 
small attraction for d < d c = d v ^ re = 2/(m — 1). This binding unbinding transition is seen 
over a characteristic length scale £y ~| v m \~ u , with v — 2/ | e m \, and e m = (m—l)(d—d^ re ). 
The length scale £y is a measure of the average separation between two contacts along the 
chain. For repulsive energy, the phase is governed by a stable fixed point where the cor- 
responding scale invariant theory is described by an infinite strength of the repulsion or 
a fermion like behavior. For d > d^ re , a threshold for the strength of the attraction is 
required for the bound state while for smaller attraction or any repulsion, the chains are 
asymptotically free. Exactly at d — d^ re , the system is at its marginal dimension where the 
repelling chains are again like free Gaussian chains but with certain log corrections, and the 
characteristic length scale diverges exponentially as the transition point is approached. 

The random system exhibits certain new features. The fact that the disorder becomes 
marginally relevant at d=l for the two body problem was observed in the context of the 
disordered wetting phenomena |L2|. The exact renormalization group (RG) analysis |5||| 



reveals, besides establishing this marginal relevance, that there is a disorder induced phase 
transition beyond the critical dimension df^f = l/(m — 1). The length scale exponent for 
this transition is found to be £f s ~ (A-Af s )~ Ur with v r = 2/(m-l)(d-d% s ), A being the 
variance of disorder. The critical value of disorder is Af s . Exactly at the critical dimension 
df^, any small disorder leads to a strong disorder behaviour over a characteristic length 
scale ~ exp(-^). However, in none of the above cases it was possible to investigate 
the strong disorder limit because of the absence of any perturbative strong coupling fixed 
point. To investigate further, a real space RG has been performed on the hierarchical lattice 



that brought out several new features [^TJ. A dynamic renormalization group approach also 
recovered these results and furthermore showed that there is no anomalous exponent for the 
free energy |7[]. This suggests that the effect of disorder is reflected in the various cumulants 
of the partition functions that characterize the phases and phase transitions [[|||. 

For the pure system, the unbound phase is characterized by the reunion partition function 
^R,m{N) which is the sum of the Boltzmann weights of all paths corresponding to the meeting 
of m chains at any point in (transverse) space at length N, the other end being tied together 
at the origin. Similarly, one can define a survival partition function Zs, m (N) for the chains 
tied at the origin at one end but located anywhere at the other end. The ratio of the two 
defines the probability of reunion. It has been shown that the reunion and the survival 
partition functions, and hence the probability, acquire anomalous scaling at the stable fixed 
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1 p^j23|| . Detailed RG analysis yields the d=l exact (vicious walker) result 
~ N~ 3 / 2 , and predicts that for marginal cases d = d^ re , the m chain reunion 

(iV) ~ N^ihiN)^ 



m 

MM 



These cases include 



point at d = 
that Z R>2 (N) 

with m body interaction has the form Z R 
the two chain problem at d = 2 and three chain at d = 1. The identical log correction, 
originating from the marginality of the interaction at the appropriate dimensions, is an 
example of "Grand Universality" for the directed polymer multicritical points This 
feature is preserved in the random case also. 

For the pure versions of cases Al, A2 and B, the unbinding critical point is at v m = 0, 
for which the chains are free. Therefore, around the transition, one can write a finite size 
scaling form 



N( m - 1 W 2 Z R , m (N)=F(x), 



(2) 



where x=\ v m \ N l l u for two chains in one dimension and x =| v m | In A" for m = d = 2 and 
m — 3, d — 1. A similar scaling form can also be written for the probabliy of reunion. A 
typical case is shown in Fig 1. 

The partition function for the chains is generated recursively in their length. For the two 
chain problem, since the randomness is only along the z direction, only the relative chain is 
considered. The center of mass is a free Gaussian chain. This transformation is not possible 
for three chains and we consider the full partition function. 

For two chains with pair interaction, the recursion relation is given by 



Z(t, t+l)= e"" [£{Z(r + e h t) + Z(t - e t , t)} + P Z(r, t)]/(2d + p), 



(3a) 



i=l 



with &i as unit vectors in the zth direction, and r denoting the position vector for the end 
point. A weight factor p has been put in for extra weightage of the relative chain. We have 
set p equal to 4. The random energy rj = v m ± A, with equal probability, whenever r is 
zero. In this relative coordinate, Z R ^(N) = Z(0,N), and Z S:2 (N) = J2 r N). For three 
chains in d = 1, 



Z(x,y,z,t + 1) 



Z(x + i,y + j,z + k,t) 

i,j,k=±l 



/8, 



(3b) 



where the Boltzmann factor for the randomness contributes only for x = y = z, i.e., only 
when three chains meet together. All chains start at origin at t — 0. For the three chain 
case, Z R ^(N) = J2 X Z(x, x, x, N), and Z Si3 (N) = J2 x , y ,z Z(x, y, z, N). 

The randomness reduces the strength of the interaction. The effective interaction v m , 
determined by 



exp(-w m ) = [exp(-v m - A) + exp(-w m + A)]/2, 



(4) 



controls the binding - unbinding transition. Since for all the three cases the pure binding- 
unbinding transition takes place at zero energy (i.e. Boltzmann factor=l), the transition 
point for a fixed v m is given by 



In cosh A r 



(5) 
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As a result a bound state can form even for v m > 0. In other words, a thermal unbinding 
is possible because randomness produces attractive pockets. The scaling of Eq. § can be 
rewritten as N (Z^ m (N)) = .F(x), where, for cases A2 and B, x = (A — A c ) In N for a given 
v m with A c satisfying Eq. [5[ 

We computed the average reunion and survival partition functions, and also the av- 
erage probability of reunion where the probability for a given realization is defined as 
P{N) = Z^ m (N)/Zs,m{N). The main distinction is that while the first two quantities may 
be considered as annealed averaging, the last one is strictly a quenched averaged quantity. 
For a given realization, the recursion relation gives exact values (upto machine precision) 
of the partitions functions. An averaging over the disorder is done by considering around 
500 samples. The lattice size is chosen such that the average transverse size of a polymer 
(~ N 1 / 2 ) does not exceed it (to avoid boundary effects). 

For a fixed strength of t>2, we plot, in Fig. 1, the data for different A, with the scaled 
variables (Z^ m (N))N as (A — A c )miV. This plot shows a nice data collapse on the pure 
curve with A c chosen as per Eq. |5| for two different values of v%. It verifies the exponential 
divergence of the length scale near the transition. We have also checked similar data collapse 
for the other cases and various values of v m . The main inference we draw from this is that 
for the average chains in presence of disorder, the nature of the transition remains the same 
as the pure system except for a shift in the critical value. This also serves as a check on the 
simulation. 

With the above knowledge, we consider a particular value A and force the system to 
be strictly at the transition point of unbinding by choosing appropriate v m . The specific 
quantity we calculate is the fluctuation in the reunion partition function or the fluctuation 
in the probability of reunion. The raw data for the cumulant (P 2 ) c have been plotted with 
TV for various A with corresponding v m determined by Eq. In case there is a disorder 
induced transition we expect the following scaling to hold good close to this critical point 

\dis 

> 

(P 2 ) c = jV-^((A - Af a )N 1/Vr ), (6) 

where, according to the RG prediction of the divergence of the length scale near the transition 
point, v r — 1 for Cases A2 and B. If we believe in Grand Universality ||, the exponent </> 
should also be the same for these two cases. The plot with (P 2 ) C N^ vs iV(A — Af s ) indeed 
shows a data collapse for a small value of Af s < .05 and <fi = 2. The two cases, A2 and B, 
for which the pure interaction is marginal, are shown in Fig 2 and 3. Numerical accuracy 
forbids simulations for very small values of A and therefore the weak disorder phase cannot 
be probed in our simulation. 

Fig 4 shows the data collapse for the case of two chains at d — 1. Since the disorder 
is expected to be always marginally relevant, the scaled variables for data collapse, from 
RG JIH, would be A log N. We obtain a data collapse for <fi = .57 ± .03. This verifies the 
exponential divergence of the length scale for the situation where the disorder is marginal. 

To conclude, we have shown that randomly interacting directed polymers undergo a 
binding-unbinding transition and, at the critical point, the fluctuation in the reunion parti- 
tion function or the probability of reunion show scaling behaviour with length scale exponents 
that agree with RG prediction. A finite size scaling form, verified for the pure case, has been 
used to study the unbinding transition.. We also obtained the exponents that describe the 
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fluctuation of the second cumulant. This exponent is the same for two chains in (2+1) 
dimensions and three chains in (1 + 1) dimensions for both of which the pure interaction is 
marginal. 

SMB and SM acknowledge the hospitality of IFF. Their visit was supported by the Indo- 
German project Number PHY-25/1. The computations were done at Julich on IBM and 
DEC alpha computers. 
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FIGURES 



FIG. 1. Data collapse for pure and random cases, two chains in 2+1 dimensions. For the pure 
case, the y— axis is NZr^{N) while for the random case it is N{Zr^ {N)). The x axis is | v 2 | log N 
for the pure case but (A — A c ) log ./V for the random case. The legend shows the values of V2 for 
the random case. 

FIG. 2. Data collapse for fluctuation of the reunion partition function. This is for two chains 
in 2+1 dimensions. The exponent = 2 and v r = 1. The inset shows the fluctuations for various 
A and length. 

FIG. 3. Data collapse for reunion probability for three chains in 1+1 dimensions. The exponents 
are = 2 and v r = 1. The inset shows the variance for various A and length. 

FIG. 4. Data collapse for fluctuation of the reunion probability for two chains in 1+1 dimen- 
sions. <p = .57 ± .03. The inset shows the variance as in Figs 2 and 3. 
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